eIF3d: A driver of noncanonical cap-dependent translation of specific mRNAs and a trigger of biological/pathological processes.

Eukaryotic initiation factor 3d (eIF3d), a known RNA-binding subunit of the eIF3 complex, is a 66 to 68-kDa protein with an RNA-binding motif and a cap-binding domain. Compared with other eIF3 subunits, eIF3d is relatively understudied. However, recent progress in studying eIF3d has revealed a number of intriguing findings on its role in maintaining eIF3 complex integrity, global protein synthesis, and in biological and pathological processes. It has also been reported that eIF3d has noncanonical functions in regulating translation of a subset of mRNAs by binding to 5'-UTRs or interacting with other proteins independent of the eIF3 complex and additional functions in regulating protein stability. The noncanonical regulation of mRNA translation or protein stability may contribute to the role of eIF3d in biological processes such as metabolic stress adaptation and in disease onset and progression including severe acute respiratory syndrome coronavirus 2 infection, tumorigenesis, and acquired immune deficiency syndrome. In this review, we critically evaluate the recent studies on these aspects of eIF3d and assess prospects in understanding the function of eIF3d in regulating protein synthesis and in biological and pathological processes.

Targeting stem-loop 1 of the SARS-CoV-2 5' UTR to suppress viral translation and Nsp1 evasion.

SARS-CoV-2 is a highly pathogenic virus that evades antiviral immunity by interfering with host protein synthesis, mRNA stability, and protein trafficking. The SARS-CoV-2 nonstructural protein 1 (Nsp1) uses its C-terminal domain to block the messenger RNA (mRNA) entry channel of the 40S ribosome to inhibit host protein synthesis. However, how SARS-CoV-2 circumvents Nsp1-mediated suppression for viral protein synthesis and if the mechanism can be targeted therapeutically remain unclear. Here, we show that N- and C-terminal domains of Nsp1 coordinate to drive a tuned ratio of viral to host translation, likely to maintain a certain level of host fitness while maximizing replication. We reveal that the stem-loop 1 (SL1) region of the SARS-CoV-2 5' untranslated region (5' UTR) is necessary and sufficient to evade Nsp1-mediated translational suppression. Targeting SL1 with locked nucleic acid antisense oligonucleotides inhibits viral translation and makes SARS-CoV-2 5' UTR vulnerable to Nsp1 suppression, hindering viral replication in vitro at a nanomolar concentration, as well as providing protection against SARS-CoV-2-induced lethality in transgenic mice expressing human ACE2. Thus, SL1 allows Nsp1 to switch infected cells from host to SARS-CoV-2 translation, presenting a therapeutic target against COVID-19 that is conserved among immune-evasive variants. This unique strategy of unleashing a virus' own virulence mechanism against itself could force a critical trade-off between drug resistance and pathogenicity.

Inhibition of SARS-CoV-2 coronavirus proliferation by designer antisense-circRNAs.

Circular RNAs (circRNAs) are noncoding RNAs that exist in all eukaryotes investigated and are derived from back-splicing of certain pre-mRNA exons. Here, we report the application of artificial circRNAs designed to act as antisense-RNAs. We systematically tested a series of antisense-circRNAs targeted to the SARS-CoV-2 genome RNA, in particular its structurally conserved 5'-untranslated region. Functional assays with both reporter transfections as well as with SARS-CoV-2 infections revealed that specific segments of the SARS-CoV-2 5'-untranslated region can be efficiently accessed by specific antisense-circRNAs, resulting in up to 90% reduction of virus proliferation in cell culture, and with a durability of at least 48 h. Presenting the antisense sequence within a circRNA clearly proved more efficient than in the corresponding linear configuration and is superior to modified antisense oligonucleotides. The activity of the antisense-circRNA is surprisingly robust towards point mutations in the target sequence. This strategy opens up novel applications for designer circRNAs and promising therapeutic strategies in molecular medicine.

Insights on the SARS-CoV-2 genome variability: the lesson learned in Brazil and its impacts on the future of pandemics.

Since the beginning of the SARS-CoV-2 spread in Brazil, few studies have been published analysing the variability of viral genome. Herein, we described the dynamic of SARS-CoV-2 strains circulating in Brazil from May to September 2020, to better understand viral changes that may affect the ongoing pandemic. Our data demonstrate that some of the mutations identified are currently observed in variants of interest and variants of concern, and emphasize the importance of studying previous periods in order to comprehend the emergence of new variants. From 720 SARS-CoV-2 genome sequences, we found few sites under positive selection pressure, such as the D614G (98.5 %) in the spike, that has replaced the old variant; the V1167F in the spike (41 %), identified in the P.2 variant that emerged from Brazil during the period of analysis; and I292T (39 %) in the N protein. There were a few alterations in the UTRs, which was expected, however, our data suggest that the emergence of new variants was not influenced by mutations in UTR regions, since it maintained its conformational structure in most analysed sequences. In phylogenetic analysis, the spread of SARS-CoV-2 from the large urban centres to the countryside during these months could be explained by the flexibilization of social isolation measures and also could be associated with possible new waves of infection. These results allow a better understanding of SARS-CoV-2 strains that have circulated in Brazil, and thus, with relevant infomation, provide the potential viral changes that may have affected and/or contributed to the current and future scenario of the COVID-19 pandemic.

Deciphering the link between Diabetes mellitus and SARS-CoV-2 infection through differential targeting of microRNAs in the human pancreas.

PURPOSE: Coronavirus Disease 2019 (COVID-19) severity and Diabetes mellitus affect each other bidirectionally. However, the cause of severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) infection on the incidence of diabetes is unclear. In the SARS-CoV-2-infected cells, host microRNAs (miRNAs) may target the native gene transcripts as well as the viral genomic and subgenomic RNAs. Here, we investigated the role of miRNAs in linking Diabetes to SARS-CoV-2 infection in the human pancreas. METHODS: Differential gene expression and disease enrichment analyses were performed on an RNA-Seq dataset of human embryonic stem cell-derived (hESC) mock-infected and SARS-CoV-2-infected pancreatic organoids to obtain the dysregulated Diabetes-associated genes. The miRNA target prediction for the Diabetes-associated gene transcripts and the SARS-CoV-2 RNAs has been made to determine the common miRNAs targeting them. Minimum Free Energy (MFE) analysis was done to identify the miRNAs, preferably targeting SARS-CoV-2 RNAs over the Diabetes-associated gene transcripts. RESULTS: The gene expression and disease enrichment analyses of the RNA-Seq data have revealed five biomarker genes, i.e., CP, SOCS3, AGT, PSMB8 and CFB that are associated with Diabetes and get significantly upregulated in the pancreas following SARS-CoV-2-infection. Four miRNAs, i.e., hsa-miR-298, hsa-miR-3925-5p, hsa-miR-4691-3p and hsa-miR-5196-5p, showed preferential targeting of the SARS-CoV-2 genome over the cell's Diabetes-associated messenger RNAs (mRNAs) in the human pancreas. CONCLUSION: Our study proposes that the differential targeting of the Diabetes-associated host genes by the miRNAs may lead to diabetic complications or new-onset Diabetes that can worsen the condition of COVID-19 patients.

Inhibition of SARS-CoV-2 Replication by a Small Interfering RNA Targeting the Leader Sequence.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has infected almost 200 million people worldwide and led to approximately 4 million deaths as of August 2021. Despite successful vaccine development, treatment options are limited. A promising strategy to specifically target viral infections is to suppress viral replication through RNA interference (RNAi). Hence, we designed eight small interfering RNAs (siRNAs) targeting the highly conserved 5'-untranslated region (5'-UTR) of SARS-CoV-2. The most promising candidate identified in initial reporter assays, termed siCoV6, targets the leader sequence of the virus, which is present in the genomic as well as in all subgenomic RNAs. In assays with infectious SARS-CoV-2, it reduced replication by two orders of magnitude and prevented the development of a cytopathic effect. Moreover, it retained its activity against the SARS-CoV-2 alpha variant and has perfect homology against all sequences of the delta variant that were analyzed by bioinformatic means. Interestingly, the siRNA was even highly active in virus replication assays with the SARS-CoV-1 family member. This work thus identified a very potent siRNA with a broad activity against various SARS-CoV viruses that represents a promising candidate for the development of new treatment options.

Higher entropy observed in SARS-CoV-2 genomes from the first COVID-19 wave in Pakistan.

BACKGROUND: We investigated the genome diversity of SARS-CoV-2 associated with the early COVID-19 period to investigate evolution of the virus in Pakistan. MATERIALS AND METHODS: We studied ninety SARS-CoV-2 strains isolated between March and October 2020. Whole genome sequences from our laboratory and available genomes were used to investigate phylogeny, genetic variantion and mutation rates of SARS-CoV-2 strains in Pakistan. Site specific entropy analysis compared mutation rates between strains isolated before and after June 2020. RESULTS: In March, strains belonging to L, S, V and GH clades were observed but by October, only L and GH strains were present. The highest diversity of clades was present in Sindh and Islamabad Capital Territory and the least in Punjab province. Initial introductions of SARS-CoV-2 GH (B.1.255, B.1) and S (A) clades were associated with overseas travelers. Additionally, GH (B.1.255, B.1, B.1.160, B.1.36), L (B, B.6, B.4), V (B.4) and S (A) clades were transmitted locally. SARS-CoV-2 genomes clustered with global strains except for ten which matched Pakistani isolates. RNA substitution rates were estimated at 5.86 x10-4. The most frequent mutations were 5' UTR 241C > T, Spike glycoprotein D614G, RNA dependent RNA polymerase (RdRp) P4715L and Orf3a Q57H. Strains up until June 2020 exhibited an overall higher mean and site-specific entropy as compared with sequences after June. Relative entropy was higher across GH as compared with GR and L clades. More sites were under selection pressure in GH strains but this was not significant for any particular site. CONCLUSIONS: The higher entropy and diversity observed in early pandemic as compared with later strains suggests increasing stability of the genomes in subsequent COVID-19 waves. This would likely lead to the selection of site-specific changes that are advantageous to the virus, as has been currently observed through the pandemic.

SARS-CoV-2 uses a multipronged strategy to impede host protein synthesis.

The coronavirus SARS-CoV-2 is the cause of the ongoing pandemic of COVID-191. Coronaviruses have developed a variety of mechanisms to repress host mRNA translation to allow the translation of viral mRNA, and concomitantly block the cellular innate immune response2,3. Although several different proteins of SARS-CoV-2 have previously been implicated in shutting off host expression4-7, a comprehensive picture of the effects of SARS-CoV-2 infection on cellular gene expression is lacking. Here we combine RNA sequencing, ribosome profiling and metabolic labelling of newly synthesized RNA to comprehensively define the mechanisms that are used by SARS-CoV-2 to shut off cellular protein synthesis. We show that infection leads to a global reduction in translation, but that viral transcripts are not preferentially translated. Instead, we find that infection leads to the accelerated degradation of cytosolic cellular mRNAs, which facilitates viral takeover of the mRNA pool in infected cells. We reveal that the translation of transcripts that are induced in response to infection (including innate immune genes) is impaired. We demonstrate this impairment is probably mediated by inhibition of nuclear mRNA export, which prevents newly transcribed cellular mRNA from accessing ribosomes. Overall, our results uncover a multipronged strategy that is used by SARS-CoV-2 to take over the translation machinery and to suppress host defences.

De novo 3D models of SARS-CoV-2 RNA elements from consensus experimental secondary structures.

The rapid spread of COVID-19 is motivating development of antivirals targeting conserved SARS-CoV-2 molecular machinery. The SARS-CoV-2 genome includes conserved RNA elements that offer potential small-molecule drug targets, but most of their 3D structures have not been experimentally characterized. Here, we provide a compilation of chemical mapping data from our and other labs, secondary structure models, and 3D model ensembles based on Rosetta's FARFAR2 algorithm for SARS-CoV-2 RNA regions including the individual stems SL1-8 in the extended 5' UTR; the reverse complement of the 5' UTR SL1-4; the frameshift stimulating element (FSE); and the extended pseudoknot, hypervariable region, and s2m of the 3' UTR. For eleven of these elements (the stems in SL1-8, reverse complement of SL1-4, FSE, s2m and 3' UTR pseudoknot), modeling convergence supports the accuracy of predicted low energy states; subsequent cryo-EM characterization of the FSE confirms modeling accuracy. To aid efforts to discover small molecule RNA binders guided by computational models, we provide a second set of similarly prepared models for RNA riboswitches that bind small molecules. Both datasets ('FARFAR2-SARS-CoV-2', https://github.com/DasLab/FARFAR2-SARS-CoV-2; and 'FARFAR2-Apo-Riboswitch', at https://github.com/DasLab/FARFAR2-Apo-Riboswitch') include up to 400 models for each RNA element, which may facilitate drug discovery approaches targeting dynamic ensembles of RNA molecules.

Analysis of SARS-CoV-2 mutations in the United States suggests presence of four substrains and novel variants.

SARS-CoV-2 has been mutating since it was first sequenced in early January 2020. Here, we analyze 45,494 complete SARS-CoV-2 geneome sequences in the world to understand their mutations. Among them, 12,754 sequences are from the United States. Our analysis suggests the presence of four substrains and eleven top mutations in the United States. These eleven top mutations belong to 3 disconnected groups. The first and second groups consisting of 5 and 8 concurrent mutations are prevailing, while the other group with three concurrent mutations gradually fades out. Moreover, we reveal that female immune systems are more active than those of males in responding to SARS-CoV-2 infections. One of the top mutations, 27964C > T-(S24L) on ORF8, has an unusually strong gender dependence. Based on the analysis of all mutations on the spike protein, we uncover that two of four SASR-CoV-2 substrains in the United States become potentially more infectious.

Genome-wide mapping of SARS-CoV-2 RNA structures identifies therapeutically-relevant elements.

SARS-CoV-2 is a betacoronavirus with a linear single-stranded, positive-sense RNA genome, whose outbreak caused the ongoing COVID-19 pandemic. The ability of coronaviruses to rapidly evolve, adapt, and cross species barriers makes the development of effective and durable therapeutic strategies a challenging and urgent need. As for other RNA viruses, genomic RNA structures are expected to play crucial roles in several steps of the coronavirus replication cycle. Despite this, only a handful of functionally-conserved coronavirus structural RNA elements have been identified to date. Here, we performed RNA structure probing to obtain single-base resolution secondary structure maps of the full SARS-CoV-2 coronavirus genome both in vitro and in living infected cells. Probing data recapitulate the previously described coronavirus RNA elements (5' UTR and s2m), and reveal new structures. Of these, ∼10.2% show significant covariation among SARS-CoV-2 and other coronaviruses, hinting at their functionally-conserved role. Secondary structure-restrained 3D modeling of these segments further allowed for the identification of putative druggable pockets. In addition, we identify a set of single-stranded segments in vivo, showing high sequence conservation, suitable for the development of antisense oligonucleotide therapeutics. Collectively, our work lays the foundation for the development of innovative RNA-targeted therapeutic strategies to fight SARS-related infections.

Global cataloguing of variations in untranslated regions of viral genome and prediction of key host RNA binding protein-microRNA interactions modulating genome stability in SARS-CoV-2.

BACKGROUND: The world is going through the critical phase of COVID-19 pandemic, caused by human coronavirus, SARS-CoV-2. Worldwide concerted effort to identify viral genomic changes across different sub-types has identified several strong changes in the coding region. However, there have not been many studies focusing on the variations in the 5' and 3' untranslated regions and their consequences. Considering the possible importance of these regions in host mediated regulation of viral RNA genome, we wanted to explore the phenomenon. METHODS: To have an idea of the global changes in 5' and 3'-UTR sequences, we downloaded 8595 complete and high-coverage SARS-CoV-2 genome sequence information from human host in FASTA format from Global Initiative on Sharing All Influenza Data (GISAID) from 15 different geographical regions. Next, we aligned them using Clustal Omega software and investigated the UTR variants. We also looked at the putative host RNA binding protein (RBP) and microRNA binding sites in these regions by 'RBPmap' and 'RNA22 v2' respectively. Expression status of selected RBPs and microRNAs were checked in lungs tissue. RESULTS: We identified 28 unique variants in SARS-CoV-2 UTR region based on a minimum variant percentage cut-off of 0.5. Along with 241C>T change the important 5'-UTR change identified was 187A>G, while 29734G>C, 29742G>A/T and 29774C>T were the most familiar variants of 3'UTR among most of the continents. Furthermore, we found that despite the variations in the UTR regions, binding of host RBP to them remains mostly unaltered, which further influenced the functioning of specific miRNAs. CONCLUSION: Our results, shows for the first time in SARS-Cov-2 infection, a possible cross-talk between host RBPs-miRNAs and viral UTR variants, which ultimately could explain the mechanism of escaping host RNA decay machinery by the virus. The knowledge might be helpful in developing anti-viral compounds in future.

Using All-Atom Potentials to Refine RNA Structure Predictions of SARS-CoV-2 Stem Loops.

A considerable amount of rapid-paced research is underway to combat the SARS-CoV-2 pandemic. In this work, we assess the 3D structure of the 5' untranslated region of its RNA, in the hopes that stable secondary structures can be targeted, interrupted, or otherwise measured. To this end, we have combined molecular dynamics simulations with previous Nuclear Magnetic Resonance measurements for stem loop 2 of SARS-CoV-1 to refine 3D structure predictions of that stem loop. We find that relatively short sampling times allow for loop rearrangement from predicted structures determined in absence of water or ions, to structures better aligned with experimental data. We then use molecular dynamics to predict the refined structure of the transcription regulatory leader sequence (TRS-L) region which includes stem loop 3, and show that arrangement of the loop around exchangeable monovalent potassium can interpret the conformational equilibrium determined by in-cell dimethyl sulfate (DMS) data.

Multiple novel non-canonically transcribed sub-genomic mRNAs produced by avian coronavirus infectious bronchitis virus.

Coronavirus sub-genomic mRNA (sgmRNA) synthesis occurs via a process of discontinuous transcription involving complementary transcription regulatory sequences (TRSs), one (TRS-L) encompassing the leader sequence of the 5' untranslated region (UTR), and the other upstream of each structural and accessory gene (TRS-B). Several coronaviruses have an ORF located between the N gene and the 3'-UTR, an area previously thought to be non-coding in the Gammacoronavirus infectious bronchitis virus (IBV) due to a lack of a canonical TRS-B. Here, we identify a non-canonical TRS-B allowing for a novel sgmRNA relating to this ORF to be produced in several strains of IBV: Beaudette, CR88, H120, D1466, Italy-02 and QX. Interestingly, the potential protein produced by this ORF is prematurely truncated in the Beaudette strain. A single nucleotide deletion was made in the Beaudette strain allowing for the generation of a recombinant IBV (rIBV) that had the potential to express a full-length protein. Assessment of this rIBV in vitro demonstrated that restoration of the full-length potential protein had no effect on viral replication. Further assessment of the Beaudette-derived RNA identified a second non-canonically transcribed sgmRNA located within gene 2. Deep sequencing analysis of allantoic fluid from Beaudette-infected embryonated eggs confirmed the presence of both the newly identified non-canonically transcribed sgmRNAs and highlighted the potential for further yet unidentified sgmRNAs. This HiSeq data, alongside the confirmation of non-canonically transcribed sgmRNAs, indicates the potential of the coronavirus genome to encode a larger repertoire of genes than has currently been identified.

Detection of bovine viral diarrhea virus genotype 1 in aerosol by a real time RT-PCR assay.

BACKGROUND: As a pestivirus of the Flaviviridae family, bovine viral diarrhea virus (BVDV), has imposed a large burden on animal husbandry worldwide, and such virus can be transmitted mainly through direct contact with other infected animals and probably via aerosols. In the present study, we aimed to develop a real-time RT-PCR method for detection of BVDV-1 in aerosol samples. METHODS: A pair of primers specific for highly conserved regions of the BVDV-1 5'-UTR was designed. The standard curve and sensitivity of the developed assay were assessed based on 10-fold serial dilutions of RNA molecular standard. The specificity of the assay was evaluated with other pestiviruses and infectious bovine viruses. The clinical performance was examined by testing 169 aerosol samples. RESULTS: The results showed that a good linear relationship existed between the standard curve and the concentration of template. The lowest detection limit was 5.2 RNA molecules per reaction. This assay was specific for detection of BVDV-1, and no amplification was found for other pestiviruses such as classical swine fever virus (CSFV), border disease virus (BDV), and common infectious bovine viruses, including BVDV-2, infectious bovine rhinotracheitis virus (IBRV), bovine parainfluenza virus type 3 (BPIV-3), bovine respiratory syncytial virus (BRSV), bovine ephemeral fever virus (BEFV) and bovine coronavirus (BcoV). The assay was highly reproducible with low variation coefficient values (CVs) for intra-assay and inter-assay. A total of 169 aerosol samples collected from six dairy herds were tested using this method. The results showed that the positive detection rate of BVDV-1 was 17.2% (29/169), which was significantly higher compared with the conventional RT-PCR. Additionally, the positive samples (n = 29) detected by real-time RT-PCR were verified by BVDV RPA-LFD, and a concordance rate of 100% was obtained between them. CONCLUSIONS: Taken together, we developed a real-time RT-PCR assay for quantitative analysis of BVDV-1 in aerosol samples, and our finding provided valuable insights into the risk on aerosol transmission of BVDV-1.